Degeneracy and Strong Fluctuation-Induced First-Order Phase Transition in the 

Dipolar Pyrochlore Antiferromagnet 
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We show that a continuous set of degenerate critical soft modes strongly enhances the first- 
order character of a fluctuation-induced first-order transition in the pyrochlore dipolar Heisenberg 
antiferromagnet. Such a degeneracy seems essential to explain the strong first-order transition 
recently observed in Gd2Sn2C>7. We present some evidence from Monte-Carlo simulations and a 
perturbative renormalization group expansion. 



I. INTRODUCTION 

In this paper we study phase transitions in Heisenberg 
magnets on the pyrochlore lattice, which consists of cor- 
ner sharing tetrahedra ' . The motivation is partly the- 
oretical, to understand the behavior of highly frustrated 
systems when there is a degeneracy, or near degeneracy, 
between different ordered states. There is also experi- 
mental motivation since experiments on Gd2Sn2 07 and 
Gd 2 Ti 2 07 have shown a rich behavior, 3,4 including multi- 
ple field transitions 5,6 which we would like to understand. 

For the family of rare earth pyrochlore systems it is 
well known 3 that dipole-dipole interactions are impor- 
tant, since the angular momentum is large (S = 7/2 for 
Gd). If one adds nearest neighbor exchange to dipole- 
dipole interactions, the Fourier transform of the total in- 
teraction J(q) is virtually independent 7 of q (and takes 
its minimum value) along the (1,1,1) directions of the re- 
ciprocal space. 3,8 This means that the magnetic ordering 
wave-vector could, potentially lie anywhere along these 
lines. 

It turns out that the phases at the endpoints, q = 
(denoted, following Ref. 10, by A) and q = (7r,7r,7r) (de- 
noted by 7r or B) are particularly important. The order- 
ing expected at q = (A-type) has been discussed by 
several authors 3,5,8 ' 10 and is shown in Fig. 1. A possible 
ordering at q = n (B-type) has been proposed in Ref. 
[10]. 

Although one would imagine that Gd2Sn2 07 and 
Gd2Ti207 should be quite similar, since the crystal struc- 
tures are the same (apart from a very small difference 
in the lattice constant), it is found that Gd2Sn2 07 has 
a strong first-order transition, 4 while Gd2Ti207 has a 
second-order transition. 12 Furthermore, Gd2Ti207 or- 
ders at q = 7T 12,13 while Gd2Sii207 appears more com- 
patible with the A phase. 14 While the small change in lat- 
tice parameter between Gd2Sn207 and Gd2Ti207 could 
change somewhat the exchange constants, it seems re- 
markable that the nature of the ordering changes so dra- 
matically. We would like to understand such a delicate 
dependence of ordering on exchange constants. 

The first order nature of the transition observed in 
Gd2Sn207 is at variance with mean- field theory which 
predicts 10 a second-order transition. In order to clarify 



FIG. 1: The A state (q = 0) of the pyrochlore lattice sta- 
bilized at low temperatures by the dipole-dipole interactions 
(Refs. [3,5,8,10]). Since the ordering is at q = 0, all tetra- 
hedra have the same spin configuration as the one shown. In 
the figure, all the spins lie onto the (xy) plane and form pairs 
of antiparallel spins that are parallel to the opposite edge of 
the tetrahedron they belong to. There are equivalent (xz) 
and (yz) states. The magnetic order is therefore character- 
ized by a n — 3-component order parameter, tj) with (xy) 
corresponding to ip = (1, 0, 0) at T=0. 



the order of the transition, and to see whether it is af- 
fected by longer range exchange interactions which lift 
the degeneracy of J(q), we have performed Monte Carlo 
simulations of the classical dipolar Heisenberg model on 
the pyrochlore lattice. These show that the apparent 
order of the transition is indeed very sensitive to the ex- 
change constants. 

We have supplemented the numerics by a perturbative 
renormalization group (RG) analysis. It is known that 
"fluctuation-induced first-order transitions" occur when 
there is no stable fixed point in a perturbative RG calcu- 
lation. This frequently occurs when the number of com- 
ponents of the order parameter n is larger than 4, 15 ~ 18 . In 
general, such transitions are expected to be only weakly 
first-order. Another microscopic scenario for a first order 
transition is the proposal of Brazovskii 19 that the exis- 
tence of a continuous set of degenerate soft modes could 
change the order of the transition. It was latter shown 
that the RG analysis of models with soft modes along spe- 
cial directions lacks stable fixed points and the models 
are indeed likely to undergo first-order transitions. 20,21 
This scenario is, for instance, relevant to the description 
of the liquid crystal transition where the degeneracy nat- 
urally comes from the isotropy of the liquid. 20 It is also 



2 



particularly relevant in frustrated magnets where frus- 
tration may, precisely, provide a large number of quasi- 
degenerate soft-modes, though with different geometri- 
cal structures. We shall consider a similar RG analysis, 
applicable for the symmetry of the dipolar pyrochlorcs 
which have lines of degeneracy, in this paper. 

In this work, we study the order of the transition of the 
classical Heisenberg model on the pyrochlore lattice with 
long-range dipolc-dipolc interactions and exchange inter- 
actions. We show, by means of Monte-Carlo simulations, 
that the transition is strongly first-order when a contin- 
uous set of soft modes is present in J(q), and becomes 
weakly first-order when that degeneracy is removed by 
including further neighbor interactions (Sec. II). In sec- 
tion III, we present a Landau-Ginzburg- Wilson simplified 
model and a perturbative renormalization group analy- 
sis of the transition, that predicts a first-order transition, 
in agreement with the numerical results in Sec. II. We 
summarize our conclusions in Sec. IV. 



We investigate ordering at q = (called A) see Fig. 1, 
and at q = tv (called B), sec Rcf. [10]. With A ordering 
the order parameter has n = 3 components and with B 
ordering it has n = 4 components as detailed in Sec. III. 
The n-component order parameters are given by 
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where is the unit vector of spin i assuming the sys- 
tem is fully ordered in component a of ordering type 
A with the (xy) state (resp. (xz), (yz)) corresponding 
to V = (1,0,0) (resp. V = (0,1,0), ip = (0,0,1)) (sec 
e.g. the arrows in Fig. 1) or B. In the course of the simu- 
lation of a finite system, the spin configuration can fluc- 
tuate between different, equivalent ordered states. We 
therefore compute the invariant quantities 



II. MONTE-CARLO SIMULATIONS 
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We present in this section some results of Monte Carlo 
simulations, using the parallel tempering approach, 22 on 
the Heisenberg model with dipolar and exchange inter- 
actions on the pyrochlore lattice. 

The classical Heisenberg Hamiltonian is: 

H = £ JijSi ■ Sj 

(hi) 



:,///;;; • ]T (^fi- - 3 ^ ' ^fr ' ^ ) 

(id) V r y r y / 
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where Si is a classical spin vector of length S = 7/2 
(for Gd 3+ ) on site i, and is the Heisenberg exchange 
between the neighbors: we will consider the first (J), 
second {J 2) and third neighbor (J 3 ) couplings. 10 

The number of spins is N — 16L 3 (L < 4) and periodic 
boundary conditions are applied. The factor of 16 arises 
because the pyrochlore lattice consists of an fee lattice 
of tetrahedra, each tetrahedron has 4 spins, and there 
are 4 sites of the fee lattice in the conventional cubic 
cell. To incorporate the B-phase with periodic boundary 
conditions, we need L to be even, so most of our results 
are for L = 2 and 4. Often the long-range dipolar in- 
teractions are cut-off beyond a couple of neighbors 23 to 
speed up the simulations. However, here we have kept a 
large number of neighbors (practically infinite) in order 
to reproduce accurately the structure of the degenerate 
states. If the dipole-dipole interaction is cut-off, ripples 
appear in the degenerate lines of soft modes 3,10 ' 11 . We 
constructed periodic repetitions of the Monte Carlo clus- 
ters and included the contributions of many blocks in 
performing the dipolar sums. Because there is no cut-off 
in the dipole interactions, the simulation becomes slow 
for large sizes, so we are limited to L < 4 (N < 1024). 



It is convenient to also compute the dimensionless Binder 
ratio 
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for both A and B ordcrings, which has the property that 
it tends to at high temperature and to 1 in an ordered 
state. (Remember that n = 3 for A-type ordering and 
n = 4 for B-typc ordering.) 

Firstly we consider the dipolar model with only nearest 
neighbor interactions. Results for the order parameters 
are shown in Fig. 2. The parameter J is taken to be the 
same as in Gd2Ti207, i.e. from the high temperature sus- 
ceptibility J = 0.4K, and the dipole-dipole interactions 
are fixed by the inter-ion distances. 3 ' 4 The order parame- 
ter associated with the A state jumps rapidly at T ~ 0.7 
K, indicating a strong first order transition at that tem- 
perature, whereas that of the B state goes to zero at 
the same temperature. These results gives a much lower 
value of the transition temperature T/v than the mean- 
field theory value of 5.3 K, though it is a bit smaller than 
the experimental value of 1 K. 

The Binder ratio shown in Fig. 3 is also strongly dis- 
continuous. In the vicinity of T c the Binder ratio gets 
negative, as expected for a first-order phase transition. 24 
Given the results of Figs. 2 and 3, we conclude that the 
transition is strongly first-order for J2 = J3 = 0. The 
very small preference for the B-phase 9-11 indicated by 
the small minimum of J(q) at q = tv, which would be 
relevant at a second order transition, is unimportant here 
because the transition is so strongly first order. 

We now include further neighbor interactions J 2 and 
J3 that lift the degeneracy of the dipolar model and select 
other states, as studied in detail in Ref. [10]. The lowest 
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FIG. 2: (Color online). Order parameter squared (xN) for 
the A phase, shown in Fig. 1, and the B phase, for J2 = J3 = 
as function of temperature T in a Monte Carlo simulation of 
the dipolar pyrochlore Heisenberg antiferromagnet. We see a 
large jump in the order parameter for the largest system size 
L = 4 signaling the onset of a strong first-order transition. 
The number of spins is given by N = 16L 3 . 
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FIG. 3: (Color online). Binder ratio for the A phase for 
different system sizes L for J2 = J3 = 0. 



part of the spectrum of J(q) is shown in Fig. 4 for q 
along the (1,1,1) direction. We see that spectrum is 
almost precisely degenerate for J2 = J3 = but that 
q = is preferred if J2 and J3 are negative, while q = 7r 
is preferred if J2 and J3 are positive. For simplicity we 
restrict ourselves to J2 = J3 and study how the character 
of the transition is modified relative to the case J2 = J3 = 
0. 

For J2 = J3 < 0, we find the same A state as for 
J-2 = J3 = 0, but the transition temperature shifts to 



J 2 =J 3 =0 
J 2 =J 3 =0.025J 
J 2 =J =-0.025J 




FIG. 4: The lowest part of the spectrum of J(q) (from Ref. 
[10]). One sees that the degeneracy of the critical soft modes 
is lifted by second and third neighbor couplings. A positive 
(resp. negative) J2 = J3 favors q = 7r (resp. q = 0). 



higher temperatures, see Figs. 5 and 6. This is expected 
from Fig. 4 since, with J2 = J3 < 0, J(q) acquires a 
well-defined minimum at q = which gets deeper with 
increasing J2 and J3. 

Even with couplings as small as J2 = J3 = — 0.061J, 
the order parameter and the Binder ratio, shown in 
Fig. 6, vary in a much more gradual way than for the de- 
generate case J2 = J3 = 0. Although only finite-size scal- 
ing on a bigger range of sizes could say whether the tran- 
sition is first or second order, it is clear that removing the 
degeneracy makes the transition less first-order compared 
with the degenerate case. These results are consistent 
with earlier simulations on model without dipole-dipole 
interactions (where larger clusters could be considered) 
which pointed out a continuous 25 or weakly first-order 26 
transition for the pyrochlore lattice with J3 < 0. 

Next we consider J 2 — J3 > which, from Fig. 4, is 
expected to favor the B-like states (q = 7r), and indeed 
this is the case as shown by Figs. 7 and 8. The mini- 
mum of J(q) at q = 7r shown in Fig. 4 is independent 
of J2 and J3 which implies that the mean-field transition 
temperature is also independent of J2 and J3. Although 
we do not have enough system sizes to attempt a serious 
estimate of Tjv using finite size scaling, it seems that the 
transition temperature is indeed quite similar to that for 
J2 = J3 = 0. 

As was also found for J2 = J3 < 0, the transition is 
much more gradual than for the degenerate case, show- 
ing that removing the degeneracy reduces the first-order 
character of the transition. The sensitive dependence of 
both the order of the transition and the nature of the 
ordered phase on J2 and J3 points to the possible rele- 
vance of these terms in explaining the difference between 
Gd 2 Sn 2 7 and Gd 2 Ti 2 7 . 

While the data for J2 — J3 — clearly indicate a 
strong first order transition, our results for J2 and J3 non- 
zero (with either sign) are not conclusive as to the order 
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order behavior found in the degenerate case. 
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FIG. 5: (Color online). Order parameters squared (xiV) for 
J2 = J3 = — 0.061J. A transition occurs to the A phase 
(q = 0) which is more gradual than for the degenerate case 
with J2 = J3 = shown in Fig. 2. 
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FIG. 6: (Color online). The Binder ratio for the A and B 
phases for different system sizes L for J2 = J3 ~ — 0.061J. 



of the transition. The smooth behavior of the Binder 
ratios shown in Figs. 6 and 8 is typical for a second order 
transition. However, we cannot rule out the possibility of 
a weak first order transition where the correlation length 
at the critical point, £ c , is large. In this case, if L < £ c the 
behavior will look like that of a second order transition. 
Only for sizes where L > £ c can one see a crossover to 
behavior expected at a first order transition. Hence for 
J2 and J3 non-zero, where the degeneracy is removed, 
the transition is either second order or weakly first order. 
However, comparing Figs. 6 and 8 with the corresponding 
figure for J 2 = J 3 = 0, Fig. 3, we see that behavior for 
J2 and J3 non-zero is very different from the strong first 
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FIG. 7: (Color online). Order parameters obtained by Monte 
Carlo simulation for J2 = J3 = 0.01J. The transition gives 
rise to the B-like state (q = 7r) and is more gradual than for 
the case of J2 = J3 = shown in Fig. 2. 
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FIG. 8: (Color online). Binder ratio for the A and B phases 
for different system sizes L for J2 ~ J3 = 0.01J. 



III. RG ANALYSIS OF THE (1,1,1) MODEL 

We now study a Landau-Ginzburg- Wilson model by 
means of the renormalization group analysis. Although 
the method is usually aimed to study second-order phase 
transitions, the lack of stable fixed points is often consid- 
ered as an indication for a first-order kind of transition. 

Given the degeneracies of the soft modes with q along 
the 4 equivalent (1,1,1) directions when J2 = J3 = (see 
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Fig. 4), the fluctuations of all these modes must be taken 
into account simultaneously. For this reason, we consider 
a model with an infinite-component order parameter (ex- 
tended to dimension d) and the fluctuations with wave- 
vectors close to these soft mode directions. The quadratic 
part of the Hamiltonian is written: 



^ ' l—l 



q- (w»-q)q 



(6) 

(7) 
(8) 



where the ii are of norm 1 and represent the i = 1, ...,4 
(1,1,1) directions given in Fig. 9. We have 8 fields, ipiq 
and (i = 1, . . . ,4), with -0 iq = V-i ,-q = V\* q - If we 
ignore the a<?jj™ term, then, when r = 0, all the modes 
with qj_,j = become simultaneously unstable. However, 
as in previous works, 20,21 we include the small dispersion 
along the (1,1,1) lines, ag 2 ™, to make the calculation well 
defined at intermediate stages. To study the degenerate 
case, it will be eliminated at the end of the calculation 
by taking m — > +oo. 



(1,-1,1) 




(1,1,1) 



FIG. 9: Lines of minimum energy in reciprocal space, given 
by the 4 equivalent (1,1,1) directions (the cube is drawn for 
convenience). The (lack of) dispersion along the (1,1,1) lines 
is shown in Fig. 4. 

The fourth-order invariants are similar to that introduced 
to describe the nematic smectic-C transitions 20 and are 
given by: 



d d q 1 d d q 2 d d q 3 d d q 4 
(27r) d (2n) d (27r) d (2n) d 



-'qi,q2,q3,q4 



H 4 (9) 



where 5 qi) q 2) q 3) q 4 ensures that the total momentum of 
the four ijj is zero. For instance, we can cancel the mo- 
mentum by choosing pairs of momenta along the same 
(1,1,1) line, i.e. by combining V*,qi with ^'-*,q2=-qu 
and similarly with q 3 and q4. That gives a fourth-order 
term Uptpitpi^i+p^i+p term, where i = 1, . . . , 4 and (i+p) 
is meant for (i + p — 4) if (i+p) > 4. Note that given 
the C4 symmetry, u p does not depend upon i, but only 
upon p = 0, 1,2. In addition, we could choose the first 
two q's along 1 and 3 for instance (see Fig. 9) and the 



other two along -2 and -4, which gives "01^2^3 ^4- An- 
other simplifying feature that we have adopted consists 
of neglecting the wave-vector dependence of the coeffi- 
cients u p . Omitting to write the qi,2,3,4 wave-vectors, 
the only fourth-order invariants are given by: 

2 4 

Hi = ^Up'^ipi'ipiipi+p'ipi+p 

p=0 j = l 

+ ^"3 (01^203^4 + V^V^) (10) 

where ^ = V-i- We call + the (1,1,1) model. 
First, we consider the Hartree correction 19 to the gap r 
(self-energy) : 



r + qi, + S 4 (r) 



(11) 



Sj(r) - ^(3u + 2 Ul +u 2 ) J ^Lg w (q) (12) 

In d = 3, if we introduce momentum cut-offs, the new 
gap r' (a — > 0) is given by: 



A 2 

r' = r + aA' ln(l + — ) 



(13) 



where a is a proportionality coefficient and A and A' the 
cut-offs. Due to the strong singularity of the right-hand- 
side, the gap does not vanish anymore. It suggests that 
the paramagnetic phase remains locally stable below the 
transition, together with other more stable phases. From 
the existence of other phases (at least at the mean-field 
level), the transition is expected to be first-order. 19 

However, the cut-offs A and A' enter explicitly the 
equation and a more controlled result can be obtained 
by the renormalization group analysis by restricting the 
integration to a shell of momentum A/6 < q < A 
with 6 > I. 20 - 21 For this we introduce new real order- 
parameters (i = 1, ...,4) : 



ipi = 4>i + i(, 



(14) 



and the quadratic and quartic terms become: 



i=l i=l 

+ u^l + 2 ) + (0 2 + 2 )][(0 2 + 2 ) + (4>\ + 2 )] 



+ u 2 [(0 2 + 2 )(0 2 _+0 2 ) + ( c 

+ "3(01020304 + 01020304) 
- H 3 (010 2 0304 + 01020304) 
+ "3(01 03 + 01 03) (^204 + 







(15) 



The derivation of the RG equations for the coupling con- 
stants is then similar to that of Ref. [20] except that we 
have to keep track of the field labels, given the anisotropy 
of Qoi in Eq. (15). By integrating over a shell of momen- 
tum A/6 < q < A, we find the recursion relations for the 
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new coupling constants: 
u' Q = 6> -(40 U 2/ 0+ 4 U 2/ + 2 U 2/ )] 
u\ = b e [ui — ( 8u 2 ii + u 3 Ii + 32uqUiIo + 8u\U 2 Iq )] 
u' 2 = b e [u2— ( + 8«2^2 + ""3^2 + 32u u 2 I a )] 
Ug = b £ [u 3 -( 8u 2 u 3 I 2 + 8uiu 3 h )] 
a' = b- 2{m -^a (16) 

where e = 4 — d and the integrals are defined by 
I * = J J^ Goi tt) goi +rW 

A. Degenerate (1,1,1) model (a — ► 0) 

The relation (16), together with the divergences of the 
integrals I p for a — > (at r = 0, finite 6) implies that a 
is a dangerous irrelevant variable for m > 1. Jo diverges 
indeed as a (d ~ 5 )/ 2 , and 7i j2 as a (d ~ 3)/2 . To take into 
account these divergences, we have to introduce rescaled 
constants vn = UicS d ~ h ^ 2 (Rcf. [20,21] and references 
therein). With a^ 5 ~ d ^ 2 Ii^ — > 0, the recursion relations 
become (we introduce X = lim a ^ a' 5 ~ d ^ 2 h)- 

u' = b m ^- d - 1 ' m ^[u -{AQul + Au\ + 2u 2 2 )I] 

u[ = b m ^- d - 1 / m ^[u 1 ~{32u Q u l +8u l u 2 )l] 

u' 2 = b m(5 - d - 1/m) [u 2 -(8u 2 1 +32u Q u 2 )2] 

u' 3 = 6^(5-^-1/™)^ (18) 

All the fixed points are unstable for d < 5 — 1/m 
(the upper critical dimension is 5 for m — ► oo) since 
u' 3 = b m ^ 5 ~ d ~ 1 ^ m ^u 3 . Although strictly speaking new 
fixed points could occur at order e 2 , the present calcula- 
tion at order e is compatible with the first-order transi- 
tion observed in the Monte-Carlo simulations. 



B. No Degeneracy 

We now remove the degeneracy, e.g. by including fur- 
ther neighbor couplings in the microscopic Hamiltonian. 

1. q = 

First of all we assume that the ordering is at q = 0. It 
is simplest to go back to the order parameter shown in 
Fig. 1, realize that there are three components ipi, i = 
1,2,3, and that the symmetry is cubic. The Hamiltonian 
of this cubic model is therefore given by 

3 

w = E t( r + q 2 )^ 2 + u o^t] + E tfti ■ ( 19 ) 

»=1 ij 

There has been a controversy regarding whether the sta- 
ble fixed point of the cubic model is the Heiscnberg or 



the cubic fixed point. A recent 6- loop expansion has 
shown that for n > 2.89, the stable fixed-point is the 
cubic one. 27 Depending on the initial values for the cou- 
pling constants, the transition could be either first-order 
or continuous. However, to stabilize the collinear states 
(with i) either (1,0,0), (0,1,0), or (0,0,1) and q = 0), 
the set of initial coupling constants leads to a first-order 
transition. 28 

We have seen that a first order transition is obtained 
both in the degenerate case and also when there is a 
well-defined minimum at q = 0. Is, then, the degener- 
acy of soft modes important or not? We note that, in 
the absence of degeneracy, the transition may be only 
weakly first-order. The problem was studied some years 
ago in the context of the pyrochlore FeF3 . For this com- 
pound, the q = state found by neutron scattering 
was characterized by a 3-component order parameter, 25 
similar to the one we have here. Monte-Carlo simu- 
lations have shown that the collinear structures with 
V> = (1,0,0); (0,1,0); (0,0,1) are preferred, but the tran- 
sition first appeared to be second-order, with unusual 
critical exponents 25 contrary to the RG argument given 
above. We can reconcile these results by suggesting that 
the transition may be weakly first-order, so that the cor- 
relation length would exceed the size of the Monte-Carlo 
cluster and the transition would appear second-order in 
the simulation. This is also comforted by a reexamina- 
tion of the Monte-Carlo results, which suggested that the 
transition is more likely to be indeed weakly first-order. 26 



2. q =77 

We now assume that the degeneracy is lifted in such a 
way that one of the four ir wave- vectors is selected. Since 
7r and — 7r are related by a reciprocal lattice vector, we 
have to take into account the fluctuations of four fields 
only, with q close to any of the 7r wave- vectors, = 
1, . . . , 4. The critical model is given by 

n = J2[(r + q 2 )^ 2 + u,^] + Ul J2 titf 

+ M3V>lV'2V'3V'4 • (20) 

This model is known to possess unstable fixed points at 
order e 2 . 18 Therefore the transition to the 7r phases is 
also expected to be first-order. 

On the basis of the LGW models alone, we would con- 
clude that the phase transitions in the dipolar pyrochlore 
are all first-order in character. Such a simple analysis 
does not say whether the transition is strongly or weakly 
first-order that is quite a relevant question when one 
comes to compare with experiments. Nonetheless, the 
results presented in this section are compatible with the 
Monte-Carlo simulations of section II. The latter are 
important, precisely to say whether the transitions are 
weakly or strongly first-order. 
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IV. CONCLUSION 

We have considered the dipolar Heisenberg model on a 
pyrochlore lattice with nearest neighbor interactions and 
a small amount of second and third neighbor interactions 
(J 2 and J3). For J 2 = J3 = the system is highly degen- 
erate, see Fig. 4, and fluctuation effects pick out ordering 
at q = (A- type). Monte Carlo simulations show that 
the transition is very strongly first order in this case, in 
contrast to mean-field theory which predicts a second or- 
der transition. A first order transition is also predicted by 
a renormalization group analysis. When the degeneracy 
is removed by including J 2 and J 3 the transition is more 
gradual, showing that the degeneracy is necessary to get 
a strong first order transition. Given the limited range of 
sizes in the Monte Carlo simulations, we cannot say from 
the simulations whether the transition is second order or 
weakly first order for J 2 = J3 0. However, according to 
a renormalization group analysis for the non-degenerate 
case, both A and B type orderings have no stable fixed 
points, indicating, presumably, a fluctuation induced first 
order transition. Usually this type of transition is only 
weakly first order, and this seems to be consistent with 
our numerical data. 

Because of the degeneracy for J 2 = J3 = 0, a small 
amount of second and third neighbor coupling can also 
change the nature of the ground state. We find that for 
Ji = <h < the A phase is retained but for J 2 = ^3 > 
we obtain a q = it (B-type) ordering. In future work 
we will study in more detail the nature of this B-type 
phase, and also consider other possible phases that occur 
when J 2 7^ Jz- It is possible that anisotropic interactions, 
suggested on the basis of high-temperature ESR 29 and by 
EPR on diluted samples, 30 ' 31 may be needed to explain 
the experimentally observed phases in detail. 



Our results provide a natural explanation for 
Gd 2 Sn 2 07 having a strong first-order transition, 4 while 
Gd 2 Ti 2 07 has a second-order transition 12 (though a 
weak first-order transition is not ruled out experimen- 
tally) ; namely second and third neighbor interactions are 
very weak in Gd 2 Sn 2 0y, but they are stronger and posi- 
tive for Gd 2 Ti 2 C>7. In this respect, ab-initio calculations 
could give some estimate of the strength of the couplings. 
This picture is also consistent with the observations that 
Gd 2 Ti 2 07 orders 12,13 at q = 7r while Gd 2 Sn 2 C"7 should 
be A-typc with equivalent sites and moments perpendic- 
ular to the local (1,1,1) directions. 14 

In the presence of a magnetic- field, 5 ' 10 Gd 2 Ti 2 07 has 
a rich phase diagram. For the future, it would also be 
interesting to perform a study of field-induced transitions 
in Gd 2 Sn 2 C>7, since this starts off with a quite different 
state in zero field. In addition to the multiple phase 
transitions expected on the basis of mean-field theory, the 
field reduces the fluctuations and, hence should reduce 
the strong first-order character of the transition. 
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